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ABSTRACT 


A Naval Postgraduate School spacecraft design class 
proposed a multiple beam, phased array, microstrip antenna as 
part of the preliminary design of a low earth orbit communica- 
tion satellite. The antenna must provide coverage over the 
satellite's entire field of view while both uplink and down- 
link operate simultaneously on the same L-band frequency. 

This thesis assesses the feasibility of the antenna 
proposed in that preliminary design. Design tradeoffs for a 
microstrip array constrained by both available surface area 
and a limited mass budget are examined. Two different 
substrate materials are considered in terms of weight and 
performance. Microstrip patch theory is applied to array 
element design and layout and antenna array theory is applied 
to determine phase and amplitude coefficients. The focus of 
the design is on obtaining the desired beam shape and orienta- 
tion, given antenna size constraints. A corporate feed method 
is discussed and a general design presented. 

Antenna performance is predicted through the use of a 
computer model based on Modal Expansion theory and results are 
plotted in a series of graphs which demonstrate the limita- 
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I. INTRODUCTION 


A. BACKGROUND 

In June 1991, Loral Cellular Systems Corporation of Palo 
Alto, California submitted a license request to the Federal 
Communications Commission for authority to build a world-wide, 
mobile communications system [Ref. 1]. Their proposed system, 
called GLOBALSTAR, is a satellite system designed to provide 
Radio-Determination Satellite Services (RDSS) for real time 
position location and tracking, as well as voice and data 
services. The system is intended to be integrated into the 
existing Public Switched Telephone Network (PSTN) and private 
cellular networks in order to provide telephone communications 
to users anywhere in the world. 

In a project sponsored by the NASA/Universities Space 
Research Association Advanced Design Program, students at the 
Naval Postgraduate School, Monterey, CA performed a prelimin- 
ary design [Ref. 2] for a satellite system based on the 
proposed GLOBALSTAR system. The focus of this thesis is on 
the design and performance analysis of a phased array L-band 
antenna, composed of microstrip patches, which was proposed as 
part of this preliminary design. The scope is limited to the 
design of the antenna panels and microstrip patches and to the 


performance analysis of overall antenna. While the beam 


forming network will be mentioned, a detailed analysis of the 
beam forming network is outside the scope of this thesis. 
Similarly, certain mechanical details such as panel hinging 
and antenna deployment mechanisms as well as specific 


fabrication methods will not be considered in detail. 


B. SYSTEM CONCEPT 

Using this system, a mobile telephone user will have 
access to world wide coverage by accessing the PSTN via 
satellite. Using a special mobile telephone unit, a user 
communicates directly with a satellite via an L-band radio 
link. The satellite receives signals from mobile users and 
retransmits these signals to an earth based gateway via an S- 
band trunk link. In the reverse process, signals from the 
gateways are received by the satellite via the S-band link and 
then retransmitted to mobile users via the L-band link. 
Gateways are distributed across the globe and are connected to 
the PSTN. 

Acquisition, synchronization, and satellite beam hand off 
coordination is provided for the gateways by a Satellite 
Operation Control Center (SOCC) as depicted in Figure 1.1. 
The information provided by the SOCC allows each gateway 
station to track and communicate with the satellites in its 
field of view.  Ephemeris data is provided to the SOCC by 
several Tracking Telemetry and Control Stations (TT&C) which 


obtain the data through satellite observations. Overall 


communications network management is provided by a Network 


- "τοι Center (NCC). 


GATE WAYS 42 


Satellite 


DATA 


TT&C GROUND STATION DISTRIBUTION 
L-Band EN (EAST) 
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C-Band ATTITUDE CONTROL 


SPACECRAFT COMMANDS 
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TT&C GROUND STATION Communications 


(WEST) νο. 
management 


Figure 1.1 System Architecture 





C. SYSTEM DESIGN 

An orbital altitude of 750 nm gives a satellite field of 
view of approximately 108 degrees (assuming a minimum 
elevation angle of 10 degrees for the satellite observer). 
This allows a constellation of 48 identical satellites in 
circular orbits to provide continuous coverage between 75 


degrees North latitude and 75 degrees South latitude. 





The constellation is divided into eight orbital planes at 
an interval of 30 degrees with six satellites in each plane. 
There is no provision for on-orbit spares. Communication 
between the satellite and gateways is accomplished via an 
uplink and downlink at center frequencies of 6533.25 MHz and 
5207.75 MHz respectively, each with a bandwidth of 16.5 MHz. 
A mobile user communicates with a satellite via an L-band link 
Operating at a center frequency of 1618.25 MHz and a bandwidth 
of 16.5 MHz. A summary of the link budget as specified by 
[Ref. 1] is contained in Appendix A. 

Each satellite provides communications coverage over its 
entire field of view through the use of six non-overlapping 
elliptically shaped spot beams. To maximize the amount of 
time a mobile user is illuminated by any one particular beam, 
beams are oriented such that their major axes are parallel to 
the spacecraft velocity vector. Figure 1.2 shows a typical 
coverage pattern for one satellite. 

In order to minimize interference between channels, the 
six beams operate in pairs such that maximum beam separation 
is achieved on the ground. This is accomplished by switching 
beams on and off in pairs so that opposite pairs are active 
for 10 ms out of a 60 ms duty cycle. During the duty cycle, 


beam pairs alternate between transmit and receive functions. 


15° FOV CONTOUR 


GROUND TRACK 





Figure 1.2. L-Band Beam Orientation 
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D. ANTENNA PERFORMANCE REQUIREMENTS 

The L-band antenna for this system must not only provide 
adequate gain, but must also correctly shape each of the six 
spot beams while maintaining low sidelobes. The antenna must 
also allow the communications payload to maintain a constant 
bit error rate (BER). Bit error rate is a function of carrier 
to noise ratio (C/N) and C/N is a function of path length. 
Because antenna beams must be statically scanned off center to 
cover the FOV, the path length for each of the six antenna 
beams is different. Path lengths vary from a minimum of 1748 
km (measured along antenna boresight) for the inner most 
beams, to a maximum of 9066 km for the outer most beams. In 
order to maintain a constant BER, an isoflux antenna is 
required. Table I-1 summarizes the L-band antenna performance 


parameters. 


TABLE I-1. ANTENNA PARAMETERS 


£,71618.25 


(MHz) 











BW (-4 dB) 25 qe 35 5 120 40 x Του 
BW (total) 209 90. 50° x 120 70 x 140 
Gain (api) | es | 55 | 40 | 


Note: Satellite FOV is 108 








II. DESIGN CONSIDERATIONS 


A. SPACECRAFT CONFIGURATION 
1. Overview 

The spacecraft is designed to be placed into orbit 
using the Delta 7925 launch vehicle. A Satellite Launch 
Dispenser (SLD) was designed so that a single launch vehicle 
can place an entire plane of satellites into orbit at one 
time. The SLD is essentially a chest of drawers where the 
spacecraft are the drawers. The spacecraft slide into the SLD 
along rails and compress a spring when properly seated. Each 
Spacecraft is then held in place by explosive bolts. At the 
appropriate time during orbit insertion, the explosive bolts 
are fired and a spacecraft is ejected into the proper orbit. 

2. Size Constraints 

The Delta payload shroud and the SLD impose serious 
size constraints on overall spacecraft size. In order to fit 
into the space provided by the SLD, the spacecraft dimensions 
can not exceed 1.4 x 1.4 x 0.6 meters. The space allowed for 
the L-band antenna alone is 1.4 x 1.4 x 0.1 meters in the 
stowed configuration. The preliminary design allowed an 
antenna mass budget of 10 kg however this proves to be 


unrealistic. 


B. ANTENNA CONFIGURATION 

In order to achieve beam steering and beam shaping 
characteristics necessary to meet the antenna performance 
requirements, a planar phased array antenna similar to the L- 
band antenna of [Ref. 3] was considered. A planar array 
allows both beam shaping and beam steering in any direction. 
For simplicity, rectangular microstrip patch array elements 
were chosen. A microstrip antenna is both lightweight and 
easily fabricated using common printed circuit board etching 
techniques [Ref. 4]. Disadvantages usually associated with 
microstrip antennas, specifically large surface area require- 
ments and narrow bandwidth [Ref. 5], were found to be insigni- 
ficant in this application. Folding array panels are employed 
SO that the antenna can be stowed in the space available and 
conveniently deployed via a spring-loaded mechanism once the 
Spacecraft is in the proper orbit. Each array panel 
corresponds to one of the six antennas necessary to form the 
appropriate beam. Figure 2.1 shows the phased array antenna 


in the deployed configuration. 


C. MICROSTRIP PATCH THEORY 
The simplest and most common type of microstrip antenna is 
the rectangular patch mounted on a dielectric substrate of 


thickness h, shown in Figure 2.2. 





Phased Array Antenna 


Figure 2.1. 


ο να. 
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SUBSTRATE FEED Fili 


WC GROUND PLANE 





Figure 2.2. Rectangular Microstrip Patch 


1. Radiation Mechanisms 
The radiation mechanism of a microstrip patch can be 
modeled in one of several ways, each with its advantages and 
disadvantages. 
a. The Wire Grid Model 
A microstrip patch can be modeled as a grid of 
fine wires [Ref. 6]. The currents on the wires are solved for 
numerically and when found, provide accurate information about 
the antenna pattern and input impedance. The disadvantage of 
this method is that it requires an extraordinary amount of 


computer resource [Ref. 7]. 


10 





b. The Modal Expansion Model 

The modal expansion model, which is valid for 
electrically thin (h«À) substrates, represents the antenna as 
a thin TM mode cavity bounded by magnetic walls [Ref. 5] 
within which the dominant TM mode is excited along with other 
non-resonant modes. As long as the dominant mode is suffi- 
ciently excited, this model provides accurate performance 
information [Ref. 7], including the effects of feed location. 
The ability to account for a feed probe and handle complex 
patch shapes are the chief advantages of this model. 

c. The Transmission Line Model 

The transmission line model considers the patch to 
be a line resonator with radiation due to fringing fields at 
the open circuited ends [Ref. 6]. The radiator consists of 
two radiating slots separated by a distance L. The slots are 
formed at the edges of the patch element, between the patch 
and the ground plane below. Field variations along the 
radiating edges are ignored. The chief disadvantage of the 
transmission line model is that it is only valid for rectan- 
gular patches. The input impedance of the transmission line 
model is also highly dependent on feed location and fails to 
accurately predict input impedance for all feed locations 
[Ref. 7]. However, formulas have been developed to predict 
input impedance for any feed location [Ref. 6]. 

Despite its drawbacks, the transmission line 


model's main advantage is, that for rectangular patches, it 


11 


provides simple design formulas which give reasonably accurate 
performance predictions. 
2. Patch Parameters 
The dimensions of a rectangular patch, as well as 
bandwidth and gain, are determined by the operating frequency 
of the antenna, the relative dielectric constant, and 
thickness of the substrate material. The following formulas 
are’ based on the transmission line model. 
a. Width and Length 
The width and length of a rectangular microstrip 


patch are given by: 














Cone 73 (25 m) 
ii MR | 
L = £ - 241 (cm) 
(2.2) 


where 


σ = speed of light (m/s), 


f, = operating frequency (MHz), 
e, = relative dielectric constant, 
e. = effective dielectric constant, 


es 2 





. 
η pee (29859 
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(£,*0.3) (240.264) 
ο A o (2.4) 


Ξ na 
(ες-0.258) (30 - 8) 


and 


h= dielectric thickness (cm) [Ref. 6]. 


b. Bandwidth 
The bandwidth of a microstrip antenna is defined 


in terms of the antenna's quality factor (Q) as follows: 


E μπε (2.5) 


Q/VSWR 
Where VSWR is less than a specified value (2:1 or 1.5:1 for 
example) and VSWR-1 at the operating frequency [Refs. 5 and 
6], the bandwidth is dependent on both the relative dielectric 


constant, € and thickness of the substrate. Thicker sub- 


p7 
Seeaces and lower values of g. give larger bandwidths. 
Typically, microstrip antennas have bandwidths on the order of 
a few percent of the operating frequency. 
c. Gain 

Gain is also affected by substrate thickness and 
relative dielectric constant. Gain is inversely proportional 
απ directly proportional to substrate thickness. To 


achieve the necessary gain, microstrip patches can be arranged 


to form linear and planar arrays. 
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D. ARRAYS 

The radiation pattern of a single array element is 
relatively wide and typically has low gain. In orders 
achieve higher directivity and additional gain, antenna 
elements can be arranged to form linear or planar arrays. 

l. Theory 

The total field of an array is given by the vector 

addition of the fields radiated by individual elements [Ref. 
B The required directivity is achieved by arranging 
individual element excitations such that fields produced by 
each element interfere constructively in the desired direction 
and interfere destructively in all other directions. The 
radiation pattern of an array is obtained by multiplying the 
field produced by a single element by an array factor (AF). 
The array factor for an array of N identical elements is given 


by [Ref. 8] as 


N 
AF = τ ο (α-ι)ψ (2.6) 
where 
N = number of elements, 
I, = element excitation coefficient, 


n 
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y = kd x cos(0) «D , (201) 


and 

kK = wave number, 

d= distance between elements (m), 

0 = direction of main beam (degrees), 

B = phase shift between elements (degrees) [Ref. 8]. 


The main beam can be scanned in any direction by 
adjusting the progressive phase shift between elements. 

M arrays of N elements can be arranged to form a 
planar array which will allow the beam to be scanned in any 
direction. Additionally, array elements can be excited 
non-uniformly to control the beam shape. If a coordinate 
system such as that of Figure 2.3 is chosen, it is shown in 


(Ref. 8] that the array factor in the X direction is given by 


M 
AF, = > Τ 187 (m-1) (kd,sin(0)cos (6$) *B,) (2.8) 


m=1 


Similarly, the array factor in the Y direction is given by 


N 
AF, x Y MC (kd,sin(0)sin($) +B,) (2 . 9) 
T 
where 
tm - excitation coefficient for m elements and 
NEL excitation coefficient for n elements. 
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Figure 2.3. Array Coordinate System 


The total array factor is the product of AL em AF, 


and is given by 


M 
TUS c y i.e) ae (kd,sin (8) cos () +B,) j (n-1) (kdysin (6) sin() +B,) 


m=1 n 


Me 


D 


n 


II 
P 


(2 hO) 

2. Beam Synthesis 
The beam in a phased array antenna is shaped by non- 
uniform excitation of the array elements. This non-uniformity 
can be in amplitude, phase, or both. Using general synthesis 
[Ref. 9] and optimization [Refs. 10 and 11] techniques, beams 


of almost any shape and contour can be formed. For less 
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stringent beam shape criteria, a simpler synthesis method 
based on the Fourier transform method presented in [Ref. 8] 
has been derived. 

In the Fourier transform method, excitation 
coefficients are related to the Fourier transform of the array 
factor expressed as a function of wy, which is itself a 


function of 0 and @. Excitation coefficients are given by 


1 π . 
- d ["AF(y) eio 
Ín E e ay 


(2.11) 
EIUS = N 
for an odd number of elements, and 
-jÍ (2m+1) jy 
In = — [arpe *? dp, 
2T J -x PENSEN 
M s m c] 
and 
-jÍ (2m-1) jy 
I, = — f AF e á αψ , 
Από -π ο 15) 
bs mM 


for an even number of elements [Ref. 8]. 


Equations 2.11 - 2.13 yield complex numbers which 
represent the amplitude and phase excitation of an element. 


The total excitation coefficient is obtained as the product of 


Ur and E 
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III. ANTENNA DESIGN 


A. MICROSTRIP PATCH DESIGN 
1. Substrate Material 

The first step in the design process is the selection 
of a dielectric substrate material. The dominant features of 
a microstrip array are controlled by substrate parameters such 
as thickness and permittivity more than by the particular 
element type [Ref. 12]. 

Because weight 45 a primary consideration in 
spacecraft design, a substrate must be chosen which not only 
has satisfactory dielectric properties, but which also has low 
density. There are a variety of substrate materials available 
with a wide range of relative dielectric constants and 
densities. Tables summarizing various materials and their 
relevant properties can be found in [Refs. 5 and 6]. 

Two representative substrate materials which were 
considered for this design are Duroid (RT 5870), produced by 
Rogers Corporation of Chandler, AZ, and Kevlar (HI 4093), 
produced by Arlon Electronic Substrates Division of Rancho 
Cucamonga, CA. 

Studies have shown [Ref. 13] that efficiency and 
bandwidth are nearly independent of patch shape  (i.e., 
circular vs rectangular) and are determined mainly by 


substrate thickness and permittivity. 
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In order to achieve good resolution in beam shape and 
steering, the number of array elements must be maximized 
within available array area and element spacing constraints. 
This is accomplished by making the microstrip patches as small 
as possible while maintaining acceptable performance charac- 
teristics. For a given operating frequency, patch size is 
determined by relative dielectric constant and substrate 
thickness. Figure 3.1 shows the effect of relative dielectric 
on patch width based on equation (2.1). Figure 3.2 shows the 
effect of dielectric thickness (h) on patch length for a given 
operating frequency and relative dielectric constant, based on 
equation (2.2). 

Duroid substrates are available with e, ranging from 
EE 10.0. Higher £. yields a smaller patch size at the 
Expense of bandwidth and gain. Rogers Duroid RT 5870 has t, 
equal to 2.23 which allows higher efficiency and bandwidth 
performance but the density of RT 5870 Duroid is 2.2 gm/cm’. 
Arlon's HI 4093 Kevlar material on the other hand, has a 
relative dielectric constant of 3.9 and a density of 1.38 
gm/ cm. 

For a given operating frequency and dielectric 
thickness, Kevlar yields a smaller patch size than Duroid and 
although efficiency and bandwidth performance are reduced, the 
results are still within acceptable parameters. Therefore, 


the lower mass of Kevlar far outweighs its disadvantages. 


1» 


Patch Width vs. Relative Dielectric 
(freq=1618.25 MHz) 


Width (cm) 


Relative Dielectric (€r) 


Figure 3.1. Effect of e. on Width of a Microstrip Patch 
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Patch Length vs. Dielectric Height 


Length (cm) 


Figure 3.2. 





πες ος 25 ΜΗΣ) 


Height (cm) 


Effect of Height and e, on Length of a 
Microstrip Patch 


ο. 


Using a solid Duroid dielectric 0.3 cm thick, the mass 
of the substrate alone would be 36.7 kg (for an antenna area 
of 1.4 m x 4.0 m). A Kevlar substrate with the same dimensions 
would yield an antenna mass of 23.3 kg. Appendix B contains 
detailed mass estimates for antenna components. 

2. Patch Dimensions 
a. Thickness 

Because the antenna substrate and ground plane 
must form the array panels, they must be made thick enough to 
provide the required rigidity. The array panels must be rigid 
enough to withstand structural loads of launch and deployment 
as well as resist warpage due to thermal cycling. A substrate 
thickness of at least 0.3 cm was selected and assumed adequate 
to meet these criteria. The antenna mass could be reduced 
without significantly affecting electrical and mechanical 
performance by supporting a thin dielectric substrate on a 
honey-comb structure such as that described in [Refs. 3 and 
4]. The effect of raising the height of the patch above the 
ground plane is to increase bandwidth [Ref. 13], while 
reducing dielectric thickness would have the opposite effect. 

b. Area 

The equations presented in Chapter II can be used 
as a starting point for determining patch dimensions which can 
later be optimized using a more sophisticated model. For an 
operating frequency of Ε-1618.25 MHz and a relative 


dielectric constant g =3.9, equations (2.1)-(2.4) yield a 


r 
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Emueh length of 5.8 cm and a width of.7.2 cm (note that the 
width is defined as the long dimension of the rectangle). The 
gain of a single patch is approximately 5 dB. Detailed cal- 
culations for length, width, efficiency, directivity, and gain 
calculations for a single patch are contained in Appendix C. 

A FORTRAN computer program which makes use of 
modal expansion theory was used in an iterative process to 
optimize the values for patch length given a constant width. 
A listing of the program is contained in Appendix F. The 
final values for length and width are summarized in Table III- 
1 along with the operating parameters for a single microstrip 


patch. 
TABLE III-1. PARAMETERS FOR A SINGLE MICROSTRIP PATCH 


€ =3. 9, h=0.3 
(cm) 
f -1618.25 MHz 


EE «τε Ie Ie 
[Modal Expansion | 4.44 | 5.9 | 5.02 ! 2.62 - 


B. PLANAR ARRAY DESIGN 





Six array panels are needed to produce the six ellipti- 
cally shaped beams. The beams are symmetrical from outboard 
to inboard with beams 1 & 6 being identical, beams 2 & 5 
identical, and beams 3 & 4 identical. The total area avail- 
able for each array panel in the deployed configuration is 


0.9338 mft. 
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1. Element Spacing 

As previously stated, the number of array elements 
needs to be maximized in order to achieve good beam 
resolution. At the same time, the array elements must be 
placed far enough apart to avoid mutual coupling which would 
degrade performance by increasing side lobe levels (SLL). It 
has been shown [Ref. 14] that for an element spacing of d > 
A/2, mutual coupling will not cause significant side lobe 
levels. Measurements made on L-band microstrip arrays have 
shown that mutual coupling levels of -25 dB or less can be 
achieved with proper element spacing [Refs. 14 and 15]. If 
the elements are separated by much more than one half 
wavelength, unwanted grating lobes are produced. An element 
spacing of about 0.55A has been shown to give good results 
[Ref. 3]. To avoid grating lobes and mutual coupling effects, 
the patch spacing for this design was chosen as 0.1017 m 
(0.55À). This spacing allows a 4 x 9 element array to be 
placed on each of the six antenna panels. 

2. Excitation Coefficients 

The Fourier transform method which was previously 
developed for a planar array was used to obtain element phase 
and amplitude excitation coefficients. Array coefficients for 
a single beam are shown in Table III-2. The resulting array 
factors are plotted in Figure 3.3 and Figure 3.4. The excita- 
tion coefficients are tabulated in Appendix D and array factor 


plots for each beam are presented in Appendix E. 
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TABLE III-2. BEAM 1 EXCITATION COEFFICIENTS 


Element T 2 3 
| Number 
d DIEN. 0.046 ` 0.046. 90.005. 
452.36 A-127.6 A-127.6 A52.36 
2 O 017 0.066. 0.066. o-oo. 
4-57.73 A129 .3 A129.3 A-57.73 
3 DD 0.083. 0.083. 0.024 
A-153.8 A28.18 A28.18 A-153.8 
ορ. 0.094. 9.094. 0.027. 
8], μὰ A-76.9 A-76.9 AWROS3. 1 
5 0.028 0. 098- 0.098 0 3928 
AO». © Aum 0. 20 A-180.0 AO.0 
"em 0.094 0.094 0.02 
-103.1 A76.9 A76.9 A-103 «2 
0.024 0.083 0.083 0.024 
A153.8 A-28.18 A-28.18 4153.8' 


Dead. 0.066 0.066 Uo 
BED. &-129.3 &-129.3' &57.73 


0.013 0.046. 0.046. ERO 3d 
A-52.36 4127.6 4127.6 &-52.36 


3. Antenna Gain 












The maximum array gain is the sum of individual 
element gains less losses associated with the array feed. The 
maximum gain is 20.6 dBi. This is well above the maximum 
required (6 dBi) for any antenna panel, however, a significant 
amount of loss will occur in the feed structure and must be 


offset by the excess gain. 
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Figure 3.3. Beam 1 Array Factor ($290) 
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Figure 3.4. Beam 1 Array Factor ($z0) 
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4. Feed Methods 

There are a variety of feed methods for microstrip 
patch arrays and for each method there exists many variations. 
In general, feed methods are grouped into three categories - 
probe, corporate, and triplate feeds - and are briefly 
discussed next. 

a. Probe Feed 

In the probe feed method, each patch is connected 

to the feed network via a coaxial-fed probe which extends 
through the ground plane and substrate. The probe is physi- 
cally attached to the patch by a solder joint. Impedance 
matching between the coax and array element can be accom- 
plished for thick substrates by designing a teardrop-shaped 
probe which cancels out probe inductance [Ref. 3] or by 
varying the diameter of the probe for thin substrates [Ref. 
16]. Obviously it would be impractical to construct an array 
where each element had its own coaxial feed cable. The weight 
and space requirements would be prohibitive for spacecraft 
applications. 

b. Corporate Feed 

A corporate feed offers a practical solution to 

the array feed problem by allowing one coaxial line to feed a 
network composed of microstrip transmission lines. A series 
of power dividers delivers power to and from the radiating 
elements. The microstrip transmission line provides equal 


path length to all elements maintaining phase between 
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elements. The two major disadvantages of the corporate feed 
are that there is not enough room on a planar array substrate 
for radiating elements and the feed network, and high line 
loss occurs due to long path lengths. The second disadvantage 
can be overcome by integrating active elements into the 
substrate [Refs. 17 and 18] or by breaking the feed network 
into subsections, each of which is fed by a low loss coaxial 
cable [Ref. 19]. 
c. Triplate Feed 

In a triplate feed (also called a hybrid feed), 
the feed network is etched on its own substrate layer and 
either laminated to the underside of the radiating element's 
substrate or attached to the backside of the ground plane. 
The former method requires the radiation of the feed network 
to be accounted for in the antenna design while the latter 
method eliminates feed radiation effects and divides the 
design problem into two parts: the array itself and the feed 
network. 

(1) Feed Network in Front of Ground Plane. A 
triplate feed located in front of the ground plane and 
directly under the radiating element's substrate simplifies 
the array construction at the expense of increased design 
complexity. The metal pins used to attach the array element 
to the feed network are eliminated through electromagnetic 
coupling. Impedance matching is controlled by the amount of 


overlap between the microstrip patch and the feed network. 
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Radiation effects can be largely reduced by locating the feed 

network as close to the ground plane as possible.  [Ref. 13] 
(2) Triplate Behind Ground Plane. If the 

corporate feed network is located behind the ground plane, 

radiation effects from the feed network are eliminated at the 

expense of a more complicated and expensive physical construc- 

tion involving plated through holes and connecting pins. 

5. Feed Design 

Although antenna feed design is discussed in [Ref. 
20], a complete design is outside the scope of this thesis. 
A brief discussion of a possible feed method is included to 
prove the feasibility of providing the required array element 
excitations. 

A triplate corporate feed located behind the antenna 
ground plane could be used. This feed method, based on those 
presented in [Refs. 6 and 19], would simplify antenna perform- 
ance analysis by eliminating feed radiation effects. Since 
the antenna is a statically scanned array, the required phase 
shifts and impedance matching could be built into the 
microstrip feed lines using microstrip transmission line 


techniques. Figure 3.5 shows one possible feed layout. 


C. SUMMARY 
An exploded view of one antenna panel is illustrated in 


Figure 3.6. 
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Figure 3.5.  Coax Corporate Feed 


PATCH SUBSTRATE 


GROUND PLANE 


vs 


| x. < di FEED SUBSTRATE 
FEED NETWORK 


Figure 3.6. Antenna Panel Assembly 
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IV. PERFORMANCE ANALYSIS 


A. ANTENNA MODEL 
Chapter II described the wire grid, transmission line, and 
modal expansion models which are used to analyze and predict 
the performance of microstrip antennas. The wire grid model 
is unsuitable for use on any but the simplest of problems due 
to the computational resource required. Because the trans- 
mission line model and the modal expansion model do not 
account for surface waves in the dielectric substrate and do 
not account for mutual coupling, they do not provide the most 
accurate solution available [Ref. 16]. However, it will be 
shown that for arrays with proper element spacing, and all but 
extreme scan angles, the modal expansion modal provides 
adequate results with relative ease of computation. The areas 
where the modal expansion model fails will be discussed and 
shown to be insignificant for this analysis. 
1. Infinite vs. Finite Arrays 
The analysis of an infinite array yields adequate 
results for the inner array elements and accounts for a change 
in impedance with varying scan angle. Infinite array analysis 
ignores edge effects due to surface waves in the dielectric. 
The analysis of a finite array includes these edge effects but 
at the expense of more difficult calculations. It has been 


shown that surface wave power increases with dielectric 
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thickness and that for a thin dielectric (hn < 0.03A), surface 
wave power is negligible [Ref. 14]. It has also been shown 
that surface wave power diminishes with increasing array size 
and that for array sizes of 7 x 7 or greater an infinite 
approximation is acceptable [Ref. 21]. 
2. Edge Effects 
Diffractions cause by the finite edge of the ground 
plane play a significant part in the radiation pattern at wide 
scan angles and in the backlobe region. The Geometric Theory 
of Diffraction (GTD) can be combined with other analysis 
methods to provide an accurate antenna performance model in 
these regions. Diffraction edge effects can be ignored for 
scan angles less than about 70' [Ref. 22]. 
3. Input Impedance 
Although the modal expansion model enables calculation 
of input impedance, patch currents near the feed point are not 
modeled. This results in inaccurate impedance. It has been 
shown that for substrate thicknesses of less than about 0.02À, 


the idealized feed model provides adequate results [Ref. 23]. 


B. SOFTWARE MODELS 
l. Method of Moments 
The method of moments solves for unknown currents on 
the surface of the microstrip patch. If a dielectric is 
present, the dielectric Green's function is used to solve the 


electric field integral equation. This method automatically 
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accounts for mutual coupling and provides accurate results for 
antenna performance. The method of moments depends on the 
accurate calculation of the elements in the impedance matrix 
for its precision [Ref. 24] and requires a large amount of 
computer resource. Several methods have been developed to 
calculate method of moment solutions more efficiently, [Refs. 
25 and 26], but comparisons with results from other models do 
not show significant differences [Ref. 24]. 
2. Modal Expansion Model 

The modal expansion model was discussed in Chapter II. 
The limitations of this model have been shown in previous 
discussions to be insignificant for this application. Because 
of its simplicity, the modal expansion model was used to 
predict the performance of this antenna. The software program 
used for this analysis was written in FORTRAN and obtained 
from the Jet Propulsion Laboratory in Pasadena, CA. A source 


code listing for this program is contained in Appendix F. 


C. MODELING SCENARIO 

Using the excitation coefficients previously obtained, 
antenna patterns for each of the six beams were generated for 
Scan angles of 0--90? to +90°, $20 and 0--90? to +90°, φ-90. 
The results were plotted in order to asses whether or not each 
beam met the design specification in terms of scan angle and 


beam width. A contour plot of the approximate ground 
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footprint of each beam was also generated in order to assess 


antenna ground coverage. 


D. RESULTS 
l. Radiation Patterns 

The radiation pattern obtained for beam 1 is plotted 
in Figure 4.1. Beam patterns for the remaining beams are 
contained in Appendix G. 

The model results demonstrate that both beam shape and 
beam orientation are attainable for beam 1. Similar results 
were obtained for beams 2 through 6 with one notable 
exception. The specification called for beams 2 and 5 to be 
scanned 30 degrees from center. When the beams were steered 
to exactly -30 degrees and 30 degrees respectively, the result 
resembled a broad side array beam pattern. Examination of the 
antenna's normal beam pattern reveals a null which prevents 
this antenna from being scanned to exactly 30 degrees. By 
adjusting the scan angle to -27 degrees for beam 2 and 27 


degrees for beam 5, acceptable beam patterns were obtained. 


34 





I 
~ 
Ul 


I 
N 
© 


I 
N 
a 


' 
UJ 
= 


ρη 
Ὅ 
G 
“τη 
(Ú 
O 
Q 
> 
“η 
3) 
© 
= 
ο 
a 


-75 -60 -45 -30 -15 0 15 30 45 60 75 90 
Theta (deg) 





Figure 4.1. Beam 1 Radiation Pattern ($-90) 


2. Pattern Separation 

Recall that the beams are required to operate in 
transmit/receive pairs such that the satellite uplink and 
downlink can be operated simultaneously at the same frequency. 
This requires adequate pattern separation to avoid mutual 
interference at the earth's surface. Figure 4.2 demonstrates 
that this design achieves adequate pattern separation for 
beams 1 and 4. Similar results were obtained for beam pair 2 


and 5 and pair 3 and 6 and are presented in Appendix C. 
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Figure 4.2. Simultaneous Beam Patterns 


3. Antenna Footprint 
In order to clearly demonstrate antenna coverage, the 
relative gain for each beam was computed over a 5 degree by 5 
degree grid centered about the satellite's nadir. The results 
for beam 2 are presented in Figure 4.3. The beam pattern is 
shown in terms of scan angle from nadir in the in-track and 
cross-track directions. Similar results were obtained for the 


remaining beams and are contained in Appendix G. 
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Figure 4.3. Beam 2 Relative dB Gain Footprint 
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V. ALTERNATIVE DESIGN 


A. ANTENNA CONFIGURATION 

An alternative design to that presented in the previous 
chapters is one which takes advantage of the fact that the L- 
band communications signal uses time division multiplexing 
(TDM). The use of TDM, as implemented in the payload for this 
spacecraft, means that each of the six antenna beams has a 60 
ms duty cycle. Since adjacent pairs of antennas are not 
operating at the same time, the antenna could be constructed 
such that there are two beams per antenna panel thus reducing 
the number of panels required. 

An integrated multiple beam microstrip array similar to 
that described in [Ref. 27] could be constructed. Two feed 
lines would be constructed on a lower substrate which would be 
electromagnetically coupled to the patch array on the 
substrate above. Each feed line would be energized at the 
appropriate time corresponding to the on-time for the 
particular beam. This configuration would simplify antenna 
construction by eliminating the need for through-hole plating 
and connecting pins between the feed network and the patch 
elements  [Ref. 28]. The deployment mechanism for this 
configuration would also be much simpler as only two hinged 
panels would have to unfold as opposed to four. There would 


also be a reduction in weight while performance would actually 
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be enhanced due to reduced losses in the smaller feed network 
required. 

Combining two antenna panels into one and using that 
single panel for two beams would almost double the number of 
patch elements available to each beam. In addition to 
increasing gain, this would allow better resolution in both 
beam shape and orientation. Increasing the number of array 
elements would also eliminate unwanted nulls which adversely 


affect scan angle. 


B. ANTENNA ORIENTATION 

The antenna in this study is oriented such that individual 
panels are perpendicular to the spacecraft's velocity vector. 
This orientation was chosen for structural and thermal control 
reasons during the early stages of the preliminary design and 
is not optimum. If the design described above were adopted, 
the antenna panels could be oriented parallel to the velocity 
vector allowing cross track beam steering to be accomplished 
by tilting two outside panels towards the center. All 
available patch elements could then be dedicated to shaping 


the beam which would result in reduced beam overlap. 
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VI. SUMMARY 


A. ANALYSIS AND DESIGN 

Two different substrate materials were considered for the 
design of a multiple beam, phased array, microstrip antenna. 
Design curves were presented to illustrate tradeoffs between 
weight and performance in the substrate material. This 
information was used to choose an appropriate substrate 
material which would provide adequate performance at a minimum 
cost in mass. 

Closed form formulae were used to arrive at a preliminary 
width and length for microstrip patch elements. A computer 
model was then used to optimize the width and length 
dimensions. Array theory was applied to derive patch 


excitation coefficients. 


B. PERFORMANCE EVALUATION 

Patch excitation coefficients were used as input to a 
computer model, which employs modal expansion theory, and 
antenna radiation patterns were generated. 

While the Method of Moments would have yielded more 
accurate results, it would have been at the expense of 
increased computer time. The Modal Expansion method used 


provided adequate results with a reasonable amount of effort. 
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The results demonstrate that the proposed design would not 
meet all the design goals. While the desired beam shapes were 
obtained, the ground footprint of the outside beams (beam 1 
and 6) show that most of the energy would be concentrated 


outside the satellite's field of view. 


C. CONCLUSIONS 

Analysis showed that this antenna could not be constructed 
within the allocated mass budget. In fact, the mass of the 
substrate material alone would be over two times the allocated 
mass budget. 

The use of rectangular microstrip patches is not an 
efficient use of the available array area. An octagonally 
shaped patch would allow more elements to be placed on a 
panel. This would allow tighter control over beam shape and 
orientation. 

The alternative design discussed in Chapter IV would 
address the short-comings in the original design and allow the 
ground footprint to be entirely within the satellite's field 
of view. However, operating two beams per antenna panel could 
require active phase shifters, increasing cost as well as 
reliability risk. 

The use of fiber optic cables instead of coaxial transmis- 
sion line and etched copper feed networks would reduce both 


weight and power loss. 
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Since the methods for determining patch dimensions 
discussed in the literature all give slightly different 
results, this antenna would have to be tested and trimmed 
after construction in order to ensure the operating parameters 


ΛΙ: 


D. FUTURE THESIS OPPORTUNITIES 

A scale model could be built and tested by increasing the 
operating frequency to yield a smaller wavelength. The 
coupled feed model discussed in the alternative design could 
also be built and tested. Various orientations, beam 
patterns, and scan angle could be explored yielding useful 


data for future phased array design projects. 
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APPENDIX A 
LINK BUDGET SUMMARY 


TABLE A-1. LINK BUDGET [REF. 1] 


E T D [ες 
Satellite to User 
— umm si πι | 0 
EMEN — 2770-31 -— 
[power toss | ao | a Power toss | 22 | o | 


Antenna Gain S/C Ant. Gain 4.0-6.0 
(isoflux) 
EN | o5 Q am [ue — (| xe | œ | 


Satellite 1389.0 km Satellite 1389.0 

Altitude- Altitude 

750 NM 

Elevation Angle degrees Elevation degrees 
Angle 


Ranse | eo | a [am | — 1389.0 


Free Space Loss Free Space 
Loss 

RX Signal -149.8 RX Signal 

Strength Strength 

Polarization Polarization 

— Loss 


[tracking toss | — oo | — db tracking toss | oo | e 


S/C Ant. Gain 4.0-6.0 Antenna Gain 
Cisof lux) 


RX Line Losses Το A Line Loss 


L-Band Amplifier L-Band 
Eff. Amplifier Eff. 


E/NoRequi red E. t / Required — 
rotat eu | — as | se [reto Ι 165 L w 
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APPENDIX B 


I. MASS ESTIMATES 


A. Kevlar Subtrate 


Density = 0.05 lb/inř Density - 1.4 gm/cm? 


The total substrate volume is given by: 


1 = 1.4 m w = 0.667 m W= 0.003 m 
Volume = l x w x h x 2 x 6 


Volume = 3.4 x 10° cw 


The total substrate mass is given by: 


MasS,,,,,. = Volume x Density 
Mass aa = 46.5 κα Μα-5,.,..-. Όψςο ο 
B. Duroid Substrate 
Specific gravity = 2.2 Density - 2.2 gm/cm? 
Total substrate mass: 
MasSy,,;,, " Volume x Density 
MasSy,,;g Ξ 74 Kg Mass, mms nter 
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e. Copper 
Jd Microstrip Patches 


Density,, = 8.96 gm/cm 


9 > 4 x 6 = 216 


Total number of elements 
Element width - 5.9 cm Element length = 4.6 cm 
Thickness,, = 050127 e 


Total volume of array elements - length x width x thickness,, 
x tot num elements 


Total volume = 74.5 cm 


Total Mass = Density x volume 


elements 


Total Mass OB κα 


elements 
2. Microstrip Feed Lines 

eres ΤΘ πη width,,., = 0.006 πι 
Volume,,, - length 


x width PS απ σος ος 


fee fee 


ae, = Density, X VOLUME eeg 
eee, = 8.9 x 10% Kg 
3. Total Copper Mass 
The total copper mass is given by: 


Massy = Mass + Mass,.a 


elements 
Mass,, - 0.701 Kg 
D. Aluminum Ground Planes 
Since a shielded triplate feed is used, two aluminum 
ground planes are needed. 
Density, - 2.7 gm/cm Sues ο πο’ 


Volume, = l xwx thickness, - 1 


Al 
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Mass, = Volume, x Density a2 


AL 


= -3 
Mass, = 1.55" > 10 Zl 


Lp. SUMMARY 

The Duroid substrate is approximately 1.5 times heavier 
than the Kevlar for a given volume. Therefore, the Kevlar 
material should be used for the substrate. The total mass for 
the antenna is estimated as follows: 


Mass - Mass 4 Massy + Mass, 


Kevlar 


A mass margin of 0.01 is chosen and the total antenna mass 


becomes 
MASS, tena = Mass x a Dind 
Mass s κα Mass δρ ~ 14 σα 
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APPENDIX C 


I. PATCH DESIGN 
Using the transmission line model, the patch parameters 


are calculated as follows: 




















Given: [Iu 618925 MHZ ε. Ξ 3.9 πο. ο απ 
l = 3(10)?/f κ - 2π/λ 
1.5 cm |o 9 
Ji 
C Es 
W = (cm) 
s W = 5.9 cm 
The effective dielectric constant is given by: 
E de 
2 Ee, = 3.6 
The patch length is then calculated as follows: 
(£,*0.3) (7+0.264) 
mel = 0.412h ND CUENCA πεν 1 σπ 
ο. ο. ο. ο. ο] 
h 
pue - 2AJ. (cm) 
Due pe Ly ee IS Ten 
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Li GAIN 
In order to calculate gain, the directivity must be 
calculated according to the following: 


AW?n? 


J. = — τν eee eee 
i [sin (60819) )2can (0) 25in (0) ae A? 


Losses, expressed as resistance, must also be found according 
to the following relations: 


Quality factor: 
i _ 3) UN 


— 


4fh 
Radiation resistance: Εξ. = 120A/W 
Copper loss: 
d T L 
ο το ο απο. —. = 2 
5 (10) προ τὰ Q 
Dielectric loss: 
ouo) no 2 
RE Q 
e, LW 
The total resistance πο Chere - R, = 379m 


With a of VSWR = 2, Bandwidth (BW) and efficiency (η) are 


found as follows: 


(VSN T 


(ο ν VSWR) 


BW - n um 
£ ας η -- 


η 


Directivity and Gain are computed as: 


5.06 dB 


D = 10 x log (Dj) D 


G 10 x log(n x Dj σΞ5δο ας 
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APPENDIX D 


EXCITATION COEFFICIENTS 


TABLE D-1. BEAM 1 EXCITATION COEFFICIENTS 


Element 

Number 
0.013 0.046 0.046 0.013 
&52.36 Dow Θ᾽ KE dig us 452.36 
0.019 On 066. 0. 066. 0.019 
A- B. 73. Aa ALS ΜΕ΄ πμ. 73. 


0.024 0.083 0.083 0.024 
&-153.8' | 428.18 &428.18' | 4-153.8' 
o 0.094. 0.094. 0.027 
4103.1 4-76.9 &-76.9 4103.1 
5 0.028 0.098. 0.098 0.028 
40.0 4180.0 4-180.0 40.0 
0920 3 0.094 0.094 ο δρυ ` 
A-103.1 476.9 476.9 A-103.1 
7 0.024. 0.083 0.083 0.024. 
| 4153.8' | 4-28.18' | 4-28.18 A153.48 | 












0.019 | 0.066. 0.066 0.019 


ONSE 0.046. 0.046. goo ` 
4-52.36 4127.6 4127.6 4-52.36 
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TABLE D-2.  BEAM 2 EXCITATION COEFFICIENTS 


| Element 
Number 
0.003 κ. 010 E 010 0.003 
4169.3 -10.68 -10.68° A169.3° 
ο. ο 0.039. ο ος O "ONES 
A-55.0 41270 412779 A-55.0 
3 0.1029" 0.099 0.099 0.029. 
&-95.3 A84.6 A84.6 A-95.3 
4 0 2045 7 ον. 02128 ο ος. 
A-137.7 A42.3 A42.3 A-137 22 


5 0.048 | 0.167 0.167 0.048. 
4-180.0 40.0 40.0 4180.0 
0.043. 0.148. 0.148. 0.043. 

4137.7 å -42.3 Á -42.3 AIRE 

0.029 0.081. 0.099. 0.029 

PE 4-84.6° &-84.6' 495.3" 





0.003. 0.010 0.010 0.003. 
Ac MES ATONGS 410.68 A- leo 
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TABLE D-3.  BEAM 3 EXCITATION COEFFICIENTS 


| Element ile 2 3 
Number 
4 0 ΌΠ 5 0.046 . 0.046 ` 0.018 
A7.7 A-172.3 A-172.3 ALS 
2 0.002 0.007. 0.007 0.002 
A5.8 A-174.2 A-174.2 A5.8 
3 Ond. s 0.083 0.083 0.024 
A-176.1 A3.8 A3.8 A-176.1 
0.051 ο. 176 9. 1946 0.051 
esL, MEM ANO 41.9 isal 
0.062 ο. 215 0. 215 0.062 
τες ο. A0.0° 40.0. ALSO. 0 
0.051 0.15 0.176 0.050 
A178.1 A-1.9 A-1.9 A]78.1 
0.024 ο. 083 ο. 083 0.024 
DC S 1... κ.α 9 CI Aad 


0.002 0.007. 0.007. 0.002 
&-5.8 4174.2. πι“ &-5.8' 


0.013 0.046. 0.046. 0.013 
ewe. | 5172 6 A23 λα 









E 


TABLE D-4. BEAM 4 EXCITATION COEFFICIENTS 


0.013 0.046. 0.046. 0.013 
4-7.7 4172% 4172.3 A-7.7 
[ e 0.002 0.007. 0.007. 0.002 
A-5.8 4174.2 4174.2 4-5.8 

3 0.024. 0.083 0.083 0.024. 
4176.1 A-3.8 4-3.8 4176.1 

4 0.051. 0.176 0.176 0.051. 
4178.1 &-1.9 A-1.9 4178.1 

5 0.062. 0.215 0.215 0.062. 
4-180.0 «0.0 40.0 4180.0 

4 0.051. 0.176 0.176 0-051, 
4-178.1 41.9 41.9 -178.1 

0.024 0.083 0.083 0.024 

4-176.1' 43.8 43.8 4-176.1 


0.002 Re 007 τ 007 0.002 
45.8 -174.2 -174.2' 45.8 


ο ου τα 046 (5 046 0.013 
A7.7 -172.3 -172.3' IE Th 
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TABLE D-5. BEAM 5 EXCITATION COEFFICIENTS 


Element 1 2 3 
Number 
0.003 0.010 0.010 κ. 003 
&-169.3' 410.68 410.68 -163 
ο. 011 0.039 a 039 omori 
A55.0' mcn “ποτ. A55.0 
3 0.029 0.099. 0.099 0.029 
495.3 A-84.6 A-84.6 495.3 
2 0-023. θ 5. 0.146 0.043 
A139. Ἢ A-42.3 A-42.3 n13777 
5 0.048 ` 0.167 0.167 0.048. 
A-180.0 A0.0 AONO A180.0 
0.043. ONES OE EM 
A-137.7 A42.3 Δ42. 3 A-137.7 
i 0.029. 0.081 0.099 0.029. 
Á-95.3 484.6 A84.6 A-95.3 


0.011. 0.039. 0.039. DRE 
&-53.0 4127.0 4127.0 &-53.0 


0.003. ο ος 0.010 | 0.003. 
4169.3 4-10.68' | 4-10.68 4169.3 
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TABLE D-6. BEAM 6 EXCITATION COEFFICIENTS 


[ΕΞ al 2 3 
Number 
HE ο ους 0.046. 0.046. 0. ΟΠ 
A-52.36 4127.6 41276 4-52.36 
2 0" 0.066 0.066 0.0m 
A457173 A-129.3 A-129.3 457.73 
3 0.024 0.083 0.083. 0.024. 
4153.8 A-28.18 A-28.18 A153.8 
3 002m. 0.094 0.094 0.027 . 
A-103.1 A76.9 A76.9 A-103.1 
5 0.028 0.098. 0.098 | 0.028 
A0.0 A180.0 A-180.0 Δ0.0 
0.027. ο ο. ο ο 0.027. 
A103.1 A-76.9 A-76.9 &103 3 
0.024 0.083 0.083 0.024 
A-153.8° A28.18° A28.18° A&-153.8 


Q0 P9 0.066. 0.066 o TOTI 
4-57.73 4129.3 4129.3 4-57.73 


0.013. 0.046 ` 0.046 | cong 
452.36 A-127.6 | Α-127.6 452.36 
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APPENDIX E 


ARRAY FACTOR PLOTS 


$008.15 230-15 
Theta (deg) 





Figure E-1. Beam 1 Array Factor ($-90) 
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Figure E-2. Beam 1 Array Factor ($20) 
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Figure E-3. Beam 2 Array Factor ($-90) 
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Figure E-4. Beam 2 Array Factor ($z0) 
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Figure E-5. Beam 3 Array Factor (4-90) 
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Figure E-6. Beam 3 Array Factor ($20) 
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Figure E-7. Beam 4 Array Factor (=90') 
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Figure E-8. Beam 4 Array Factor ($20) 
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Figure E-9. Beam 5 Array Factor (ΦΞ90) 
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Figure E-10. Beam 5 Array Factor ($z0) 
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Figure E-11. Beam 6 Array Factor ($-90) 
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Figure E-12. Beam 6 Array Factor ($-0) 
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APPENDIX F 


DESIGN.FOR 


.MICROSTRIP ANTENNA DESIGN CHARACTERISTICS CALCULATED BY 


MODAL 


. THEORY, ALL INPUT LENGTH ARE IN CENTIMETER, BY DR. HUANG 


COMPLEX CJ,YWY,AY,DEL(5),C1,C2,D4,SKY,FC,WC 
REAL A,B,EPS 

CJ=(0.,1.) 

PI-3.14159265 

TPI-2.*PI 

DPR-180./PI 


TOL=.0001 


WRITE (6,11) 

FORMAT (' CENTER FREQ=?') 

READ (5,99) FREQ 

WRITE (6,1) 

FORMAT (' WIDTH OF PATCH=?') 

READ (5,99)A 

FORMAT (F5.8) 

WRITE(6,2) 

FORMAT(' RESONANT LENGTH OF PATCH=?' ) 
READ(5,99)B 

WRITE (6,3) 

FORMAT (' DIELECTRIC CONSTANT=? ') 

READ (5,99) EPS 

WRITE (6,4) 

FORMAT (' DIELECTRIC THICKNESS=? ') 
READ (5,99)T 

WRITE (6,5) 

FORMAT (' PROBE LOCATION=?') 

READ (5,99) YO 

WRITE (6,6) 

FORMAT (' LOSS TANGENT=? ') 
READ (5, 99) TLOS 

FF=1./SQRT(1.+10.*T/A) 
EPE=(EPS+1.)/2.+(EPS-1.)/2.*FF 

FF= (A/T+0.262) / (A/T+0.813) * (EPE+0.3) / (EPE-0.258) *0.412 
WAMDA=2 .1*B*SORT (EPS) 
GW20.00836*A/WAMDA 
BW20.01668*FF*A*EPE/WAMDA 
FY20.774740.5977* (A/B-1.) -0.1638* (A/B-1.) **2. 
YWY=CMPLX (GW, BW) 
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20 


41 


42 


a1 


a2 


332 


34 


37 


AY-CJ*TPI*376.7/WAMDA*T/A*YWY*FY*B 
DEG DS (0 TOR 

DO 20 I-1,4 
C122.*AY* (CMPLX (PI,O0.) -DEL(I)) 


C2=AY*AY+2 .*DEL (I) *PI-DEL (I) *DEL (I) -CMPLX (PI*PI,0.) 


DEL (I+1) =C1/C2- (DEL (I) ) **3./3. 
CONTINUE 

D4=DEL (5) 

SKY= (CMPLX (PI, 0.) -D4) /B 

FC=CMPLX (1.,0.) -CJ*TLOS 
WC23.E10/SQRT (EPS*FC) *SKY 

FR=REAL (WC) /TPI 
OR-REAL (WC) /2. /AIMAG (WC) 

FF=COS (PI*YO/B) 

CC=EPE*8 .854*1.E-14*A*B/2./T/FF/FF 
RAD-QR/REAL (WC) /CC 
RCU=0.00027*SORT (FR/1.E+9) *B/A*QR*OR 
RDI-230.*TLOS*T*WAMDA/EPS/A/B*QR*QR 
RIN=RAD+RCU+RDI 

Q=RIN*REAL (WC) *CC 

DFR=FR/Q 

EFF=RAD/RIN 

FR=FR/1.E+9 


DFREQ-FREQ-FR 
IF((ABS(DFREQ)).LE.TOL)GO TO 42 
WRITE (6,41) FR, DFREO 
IF(FR.LT.FREQ) THEN 

B-B-.0001 
ELSE 

B=B+.0001 
END IF 
GOTOO 


FORMAT (' FREQUENCY=',F10.5, 3X, 'DEPREOS E OR 


WRITES 2 DER 


FORMAT (' RESONANT FREQUENCY IN GIGAHERTZ-',F10.5) 


WRITE (6,32) RIN 

FORMAT (' INPUT RESISTANCE-',F10.5) 
DFR=DFR/1.E+6/2. 

WRITE (6,33)DFR 


FORMAT (' BANDWIDTH IN MEGAHERTZ=',F10.5) 


WRITE (6,34) EFF 


FORMAT (' RADIATION EFFICIENCY=',F10.5) 


PERCBW=1/ (Q*SQRT (2.0) ) 
WRITE (6,37) PERCBW 
FORMAT (' % BW=',F10.5) 
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38 


39 


OQ OQ araa 


CM 
€. 


Ω 


29 
98 


WRITE (6,38)RCU,RDI,RAD,Q 

FORMAT ('RCU=',F10.5,2X, 'RDI=',F10.5,2X, 'RAD=',F10.5, 
2X, 'O=', F10.5) 

WRITE (6,39)A,B 

FORMAT ('WIDTH=! , F10.5,5X, 'LENGTH=' , F10.5) 


STOP 
END 


MICARY.FOR 


..LINEARLY POLARIZED PLANAR ARRAY FORMED BY LINEARLY 

. . POLARIZED MICROSTRIP ELEMENTS 

..WITH CROSS-POL, MODAL EXPANSION MULTIMODES ARE USED. 
..NO GTD DIFFRACTION, ALL DIMENSIONS IN WAVELENGTHS 

. PROGRAMMED BY DR. J. HUANG AT JET PROPULSION LAB. 


DIMENSION DBR (362) ,DBL(182) , XAX (362) ,XX(2),YY(2), 
LIANG (40) 

DIMENSION PHASR(182),PHASL(182),XE(40), YE(40), 
1AMP (40), PHAS(40) 

COMPLEX CJ,CJJ,EX,EY,EZ,X0,Y0O,2Z0,ETH, EPH, ER, EL 
BSUPLEX FX,FY,FZ,X1,Y1,Z21,X2,Y1Y2,22,X3,Y3,23 
COMMON/DD1/EPS,T,A,B,I2,IP,XS,YS,ZS,CJJ 
COMMON/DD2/CJ,PI,TPI,DPR 

COMMON/DDA /THPR , PHPR, RR 


MODIFIED TO READ DATA FROM A GEOMETRY FILE 
BY BILL BARFIELD, NPGS 

OPEN (33, FILE='DATA') 

OPEN (35, FILE='GEOMETRY ' ) 


ο ο. 1.) 

IJ (].E-15,1.E-15) 

PI=3 .14159265 

TPI=2.*PI 

DPR-180./PI 

WRITE (6,1) 

FORMAT (' DIELECTRIC CONSTANT AND THICKNESS=? ' ) 


READ (5,99) EPS,T 
READ (33, *)EPS,T 
ΤτΕ, 99) ΕΡΤ 


FORMAT( F8.6,F8.6) 

FORMAT (I8) 

WRITE(6,2) 

FORMAT(' MICROSTRIP WIDTH AND LENGTH-?!) 
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τος 


READ(5,99)A,B 
READ (33,*)A,B 
WRITE (6,99)A,B 


..A,B ARE THE PHYSICAL WIDTHS OF THE METALLIC PATCH 


..A ALONE X-AXIS, B ALONE Y-AXIS, X-AXIS IS HORIZONTAL 
WRITE (6, 3) 
FORMAT (' FIELD POINT DISTANCE AND PATTERN CUT 
/ANGLE=? ' ) 


READ (5,99) RF, PHD 
READ (33, *) RF, PHD 
WRITE (6,99) RF, PHD 
WRITE (6, *) RF, PHD 


WRITE (672) 
FORMAT (' FEED PROBE DIAMETER & DISTANCE FROM EDGE=?' ) 


READ (5,99)D,YP 
READ (33, *) D, YP 
WRITE (6,99) Dp we 


WRITE (6, 4) 
FORMAT (' NUMBER OF ELEMENTS=?' ) 


READ (5,98)N 
READ (35,*)N 
WRITE (6,*)N 


WRITE (6,6)I1 
FORMAT (' NUMBER! ,1I2,2X,'X,Y LOCATION, AMPLITUDE, PHASE & 
/ORIENTATION=? ' ) 


DO του ΤΞΕ Ἢ 


READ (5, *) XW, YW, AMG, DEG, IANG(I) 
READ (35, *) XW, YW, AMG, DEG, IANG (I) 
WRITE (6, *) XW, YW, AMG, DEG, IANG (I) 


XE (I) =XW 

YE yw 

AMP (I) =AMG 

PHAS (I) =DEG/360.*TPI 
CONTINUE 

Dz5.*D 
SK2=EPS*TPI*TPI 
PIB=(PI/B) **2 
PIA- (PI/A)**2 

Μαη AVES 
Fl=1.E-7 

CYP=COS (PI*YP/B) 
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F2=2./A/B/ (SK2-PIB) *CYP 
SD=SIN(PI*D/A) / (PI*D/A) 

F3=4./A/B/ (SK2-4.*PIA-PIB) *CYP*SD 
CYP=COS (TPI*YP/B) 

F4=4./A/B/ (SK2-4.*PIA-4.*PIB) *CYP*SD 
F5-2./A/B/ (SK2-4.*PIB)*CYP 

F6=2./A/B/ (SK2-4.*PIA) *SD 

F7=2./A/B/ (SK2-9.*PIB) *COS (3.*PI*YP/B) 


C PRITE(6,7)F1,F2,F3,F4 
7 ο ewer Oe S,5x,'RF2=',F10.5,5X, 'F3=',F10.5, 5X, ' 
HET-' F10.5) 
€ WRITE(6,8)F5,F6,F7 
8 EE Wu NES UON 5X '"R6z',;RFI0.575X,'F72',F10.5) 


De 200 I=1,181 

THDD=FLOAT (I-1) /2.-45. 

THD=ABS (THDD) 

THR=THD/DPR 

PHR=PHD/DPR 

IF (THDD.LT.0O.) PHR=PHR+PI 

XS=RF*SIN (THR) *COS (PHR) 

YS=RF*SIN (THR) *SIN(PHR) 

ZS=RF*COS (THR) 

FX-CJJ 

FY-CJJ 

FZ-CJJ 

DO 104 I4=1,N 

EX-CJJ 

EY-CJJ 

EZ-CJJ 

DO 102 I2-1,4 

XA- (B4T) /2. 

IF(I2.EQ.2)XA--XA 

IF(I2.EQ.3) XA- (A4T) /2. 

IF(I2.EQ.4) XA-- (A4T) /2. 
-XS-XE(I4)-XA 

Y=YS-YE (I4) 

Z-ZS 

IF(I2.LT.3)X-2XS-XE(I4) 

IF(I2.LT.3)Y-YS-YE(IA4)-XA 

RR-SQRT (X*X4Y*Y4Z*Z) 

XU-X/RR 

YU=Y/RR 

ZU=Z/RR 

IF (ABS (XU) .LT.1.E-8) XU=1.E-8 

IF (ABS (YU) .LT.1.E-8) YU=1.E-8 

| PHPR=ATAN2 (YU, XU) 

IF (PHPR.LT.0O.) PHPR=PHPR+TPI 

IF(I2.LT.3)PHPR-PHPR-PI/2. 

IF (PHPR.LT.0O.) PHPR=PHPR+TPI 

THPR=ACOS (ZU) 
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41 


42 


43 


50 


5 1 


52 


53 


60 


61 


CALL SLOT4 (X0,Y0,Z0,X1 Y1 ρι ρου. 
IF(IANG(IA4).EQ.270)GO TO 50 
IF(IANG(I4).EQ.90)GO TO 60 
IF(IANG(I4).EQ.180)GO TO 70 

IF(I2.EO.2) G0 URS Nd 

IF(I2.EQ.3)GO TO 42 

IF(I2.EQ.4)GO TO 43 

EX=EX- (F1-F2+F5-F7) *YO- (F4-F3+F6) *Y2 

EY=EY+ (F1-F2+F5-F7) *X0+ (F4-F3+F6) *X2 

EZ=EZ+ (F1L-F2+F5-F7) *Z0+ (F4-F3+F6) *Z2 

GO TO 102 

EX=EX+ (F2+F1+F5+F7) *Y0+ (F3+F4+F6) *Y2 

BY=BY- (F2+F1+F5+F7) *X0- (F3+F4+F6) *X2 

BZ=EBZ- (F2+F1+F5+F7) *Z0- (F3+F4+F6) *Z2 

GO TO TO? 

EX=EX+ (F1-F6) *X0+ (F2-F3) *X1+ (F4-F5) *X2-F7*X3 
EY-EY-4 (F1-F6) *Y0+ (F2-F3) *Y1+ (F4-F5) *Y2-F7*Y3 
EZ=EZ+ (F1-F6) *Z0+ (F2-F3) *Z1+ (F4-FS) *Z2-F7*Z3 
GO TO το) 

EX=EX- (F1-F6) *X0+ (F3-F2) *X1- (F4-F5) *X2+F7*X3 
EY=EY- (F1-F6) *Y0+ (F3-F2) *Y1- (F4-F5) *Y2+F7*Y3 
EZ=EZ- (F1-F6) *Z0+ (F3-F2) *Z1- (F4-F5) *Z24+F7*Z3 
GO TO 102 

IF(I2.EQ.2)GO TO 51 

IF (I2.EQ.3)GO TO 52 

IF(I2.EQ.4)GO TO 53 

EX=EX- ( (F1-F6) *Y0+ (F3-F2) *Y1+ (F4-F5) *Y2+F7*Y3) 
EY=EY+ ( (F1-F6) *X0+ (F3-F2) *X1+(F4-F5) *X2+F7*X3) 
EZ=EZ+ ( (F1-F6) *Z0+ (F3-F2) *Z1+ (F4-F5) *Z2+F7*Z3) 
GO TO 102 

EX=EX- (- (F1-F6) *Y0+ (F2-F3) *Y1- (F4-F5) *Y2-F7*Y3) 
EY=EY+ (- (F1-F6) *X0+ (F2-F3) *X1- (F4-F5) *X2-F7*X3) 
EZ=EZ+ (- (F1-F6) *Z0+ (F2-F3) *Z1- (F4-F5) *Z2-F7*Z3) 
σο το 102 

EX=EX+ ( (F1-F24+F5-F7) *X0+ (F4-F3+F6) *X2) 

EY=EY+ ( (F1-F2+F5-F7) *Y0+ (F4-F3+F6) *Y2) 

EZ=EZ+ ( (F1-F2+F5-F7) *Z0+ (F4-F3+F6) *Z2) 

(CONSO IN 

EX=EX+ ( (-F2-F1-F5-F7) *X0+ (-F3-F4-F6) *X2) 
EY=EY+ ( (-F2-F1-F5-F7) *Y0+ (-F3-F4-F6) *Y2) 
EZ=EZ+ ( (-F2-F1-F5-F7) *Z0+ (-F3-F4-F6) *Z2) 

ο πο ο 

IF(I2.EQ.2)GO TO 61 

IF(I2.EQ.3)GO TO 62 

IF(I2.EQ.4)GO TO 63 

EXZEX- ((Fl-F6)*YO-«(F2-F3)*Y1-«(FA4-F5)*Y2-F7*Y3) 
EY=EY+ (F1L-F6) *X0+ (F2-F3) *X1+ (F4-F5) *X2-F7*X3 
EZ=EZ+ (F1-F6) *Z0+ (F2-F3) *Z1+ (F4-F5) *Z2-F7*Z3 
COLTO 102 

EX=EX+ (F1-F6)*Y0+ (F2-F3)*Y1l+ (F4-F5)*Y2-F7*Y3 
EY=EY- ( (F1-F6) *X0+ (F2-F3) *X1+ (F4-F5) *X2-F7*X3) 
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EZ-EZ-((Fl-F6)*Z0O-«(F2-F3)*Z1-«(F4-F5)*Z2-F7*Z3) 
ESTO 102 
62 BEX=EX+ (F1+F2+FS5+F7) *X0+ (F3+F4+F6) *X2 
BY=EY+ (F1+F24+F5+F7) *Y0+ (F3+F4+F6) *Y2 
EZ=EZ+ (F1+F2+F5+F7) *Z0+ (F3+F4+F6) *Z2 
GO TO 102 
63 EX=EX+ (F2-F1-F5+F7) *X0+ (F3-F4-F6) *X2 
EY=EY+ (F2-F1-F5+F7) *Y0+ (F3-F4-F6) *Y2 
EZ=EZ+ (F2-F1-F5+F7) *Z0+ (F3-F4-F6) *Z2 
ESTO 102 
70 NEI2.EO0.2)GO TO 71 
NET2.EQ.3)GO TO 72 
IF(I2.EQ.4)GO TO 73 
EX=EX- ( (F1+F2+F5+F7) *Y0+ (F3+F4+F6) *Y2) 
BY=EY+ (F1+F2+F5+F7) *X0+ (F3+F4+F6) *X2 
EZ=EZ+ (F1+F2+F5+F7) *Z0+ (F3+F4+F6) *Z2 
GO TO 102 
71 EX=EX- ( (F2-F1-F54+F7) *Y0+ (F3-F4-F6) *Y2) 
EY=FY+ (F2-F1-F5+F7) *X0+ (F3-F4-F6) *X2 
EZ=EZ+ (F2-F1-F54+F7) *Z0+ (F3-F4-F6) *Z2 
@® TO 102 
72 EX=EX+ (F1-F6)*X0+ (F3-F2)*X1+(F4-F5)*X2+F7*X3 
EY=EY+ (F1-F6)*Y0+ (F3-F2)*Y1+(F4-F5)*Y2+F7*Y3 
EZ=EZ+ (F1-F6) *Z0+ (F3-F2) *Z1+ (F4-FS) *Z2+F7*Z3 
ee, TO 102 
73 EX=EX+ (F6-F1) *X0+ (F2-F3) *X1+ (F5-F4) *X2-F7*xX3 
EY-EY-« (F6-F1)*YO-(F2-F3)*Y1-(F5-F4)*Y2-F7*Y3 
EZ-EZ-(F6-F1)*ZO-(F2-F3)*Z1-(F5-F4)*Z2-F7*Z3 
102 CONTINUE 
FX-FX4EX*AMP (I4) *CEXP(-CJ*PHAS (I4) ) 
FY=FY+EY*AMP (14) *CEXP (-CJ*PHAS (I4) ) 
FZ=FZ+EZ*AMP (14) *CEXP (-CJ*PHAS (I4) ) 
104 CONTINUE 
ETH=FX*COS (PHR) *COS (THR) +FY*SIN (PHR) *COS (THR) -FZ*SIN 
/ (THR) 
EPH=-FX*SIN (PHR) +FY*COS (PHR) 
ο ER- (ETH4«CJ*EPH) / SQRT (2.) 
G EL= (ETH-CJ*EPH) /SQRT (2.) 
ER-ETH 
EL-EPH 
AER-CABS (ER) 
AEL-CABS (EL) 
IF(AER.LT.1.E-15)AER-1.E-15 
IF(AEL.LT.1.E-15)AEL-1.E-15 
IF (AER.LT.1.E-8) ER=CMPLX(1.,1.) 
IF (AEL.LT.1.E-8) EL=CMPLX(1.,1.) 
DBR (I) =20.*ALOG10 (AER) 
DBL (I) =20.*ALOG1O (AEL) 
PHASR (I) =ATAN2 (AIMAG (ER) , REAL (ER) ) *DPR 
PHASL (I) =ATAN2 (AIMAG (EL) , REAL (EL) ) *DPR 
XAX (I) =THDD 
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WRITE (6, *) XAX, DBR, DBL 


CONTINUE 
YMAXL=-1000. 
YMAXR=-1000. 
DO 22 I=1,181 
IF (YMAXR.LT.DBR (I) ) YMAXR=DBR (I) 
IF (YMAXL.LT.DBL (I) ) YMAXL=DBL (I) 
CONTINUE 
YMAX=YMAXR 
IF(YMAXR.LT.YMAXL)YMAX-YMAXL 
DO 21 I-1,181 
DBR (I)-DBR(I)-YMAX-2. 
DBL (I) =DBL (I) -YMAX-2. 
IF (DBR(I) .LT.-40.)DBR(I)=-40. 
IF (DBL(I) .LT.-40.)DBL(I)=-40. 
CONTINUE 
CABLE BGNEEE 
CALL PLFORM('LINLIN!,7.43,4.6) 
XX (1) =-180. 
XX (2) 2180. 
YY (1)2-40. 
YY (2) —0- 
CALL PLSCAL (XX, 2,080908, YY, 2, 040504) 
CALL PLGRAF('LP MICROSTRIP ARRAY', 'ANGLE', 'DB') 
CALL PLCURV (XAX,DBR,181, 0,0) 
CALL PLNUP 
CALL PLNTYP (5) 
CALL PLNDN(0.,0.) 
CALL PLCURV (XAX,DBL,181, 0,0) 
CALL ENDPLT 


J=181 

T=-181 

DO 211 1-15 362 

XAX (I)z(I-1)/2.-45. 
IF(XAX(I).GT.90.) XAX (I) -XAX (I) -180. 
DBR (I) =DBR (J) 

J=J-1 

CONTINUE 


OPEN (38, FILE='RESULTS.PRN' ) 


WRITE (6,699) 
FORMAT (' WRITING RESULTS TO FILE...') 


FORMAT (F8.3,1X,F8.3,1X) 


Do σοι τι, 5ο 
WRITE (38, *) XAX(I) , DBR(T) 
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701 


1000 


CONTINUE 
WRITE(6,*)' DONE...' 
CLOSE (33) 
CLOSE (35) 
CLOSE (38) 


STOP 
END 


PesROULINE SLOT4 (EX, EY, EZ, EX1, EY1, EZ1, EX2, 


Ww EZ2,EX3,EY3,E23) 


C...RADIATION FROM DIFFERENT MODE SLOTS 


HT 
T2 





COMPLEX CJ,FAC,EPH,ETH,EX,EY,EZ,CJJ 
COMPLEX EX1,EY1,EZ1,EX2,EY2,EZ2,EX3,EY3,E23 
COMMON /DD1/EPS,T,AA,BB,12,IP,XS,YS, ZS, CJJ 
COMMON /DD2/CJ,PI,TPI,DPR 

COMMON /DD4/THR, PHR,R 

FAC=CEXP (-CJ*TPI*R) /R 


A=T/2. 

B= (AA+T) /2. 

IF(I2.GT.2) B=(BB+T) /2. 
CPHZCOS (PHR) 

SPH=SIN (PHR) 

CTH=COS (THR) 


STH=SIN (THR) 

ARG1=TPI*A*CPH*STH*SQRT (EPS) 
ARG2-TPI*B*SPH*STH 
NUUABS(ARGI).LT.1.E-4)GO TO 11 
F1-SIN(ARG1)/ARG1 

€O TO 12 

Fl=1. 

IF(ABS(ARG2).LT.1.E-4)GO TO 13 
F2-SIN(ARG2) /ARG2 

COTTO 14 

F2-1. 

JP (ABS (ABS (ARG2) -PI/2.) .LT.1.E-4)GO To 15 
F3=ARG2*COS (ARG2) / (ARG2*ARG2-PI*PI/4.) 

πο το 16 

F3=-1./PI*ARG2 

F4-2COS (ARG2) / (ARG2*ARG2-PI*PI) 
NUUABSUABS(ARG2)-1.5*PI).LT.1.E-4)GO TO 17 
F5-ARG2*COS (ARG2) / (ARG2*ARG2-PI*PI*9./4.) 
GO TO 18 

F5--1./2. 

EPH--CJ*F1*F2*SPH*CTH*FAC 
ETH-CJ*F1*F2*CPH*FAC 
EX-ETH*CTH*CPH-EPH*SPH 
EY-ETH*CTH*SPH4EPH*CPH 

EZ--ETH*STH 
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ΒΡΗΞΕ-ΡΒΙ ΕΙΠΕ πα. 
BTH=F1*F3 *CPH* FAC 
EXI-ETH*CTH*CPH-EPH*SPH 
HY 1=ETH* CTH* SPH+EPH* CPH 
EZ1=-ETH*STH 
EPHZCO*PI*Fl*FA*SPH*CTH*FAC 
ETHZ-CO*PEITBEICUHBACCDBLBAS 
EX2=ETH* CTH* CPH-EPH*SPH 
EY 2=ETH* CTH* SPH+EPH* CPH 
EZ22-ETH*STH 
EPH=F1*F5*SPH*CTH*FAC 
ETH=-F] σε s Dp G 
EX3ZETH*CTH*CPH-EPH*SPH 
EY3SEIH*CIH*SPHESEPH*CBH 
EZ3--ETH*STH 
RETURN 

END 


III.  GEOM.FOR 


C...ROUTINE TO CREATE GEOMETRY FOR MICARY.FOR 
DIMENSION X(4),Y(9),SMAGD(9,4) , SPHASED(9, 4) 
REAL L,W,AMP,PHASE,LAMBDA,FREQ, SPACE 


OPEN (5, FILE='GEOMETRY ' ) 
OPEN (33, FILE='MAG' ) 
OPEN (35, FILE=' PHASE' ) 


C...READ IN MAGNITUDE AND PHASE DATA FROM FILE 


NUMX=4 
NUMY=9 
NUMTOT=NUMX*NUMY 


DO 4 I=1,NUMY 

READ (33, *) SMAGD(I, 1) , SMAGD(I, 2) , SMAGD(I,3) , SMAGD(I, 4) 
READ (35, *) SPHASED(I,1),SPHASED(I,2) ,SPHASED(I,3), 
/SPHASED (I,4) 

WRITE (6, *)SMAGD(I,1),SMAGD(I,2),SMAGD(I,23), SMAGD(I, 4) 
WRITE(6,*)SPHASED(I,1),SPHASED(I,2),SPHASED(I,3), 
/SPHASED (I,4) 


3 CONTINUE 
12 FORMAT (F10.5,F10.5 PTOS ETON 
C...FREQUENCY IN HZ 


FREQ=1618 .25E6 
LAMBDA=3 .E8/FREQ 
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C...DIMENSIONS IN CM 
L-4.44 
Wz5.9 
AMP=1.0 
PHASE=. 0 
ORIENT=0. 
SPACE=.55*LAMBDA 


C...CONVERT TO WAVELENGTHS 
L-L/100. 
L-L/LAMBDA 
W-W/100. 
W=W/LAMBDA 


DO 10 I=1,NUMX 
X (I) =(I-1) *SPACE+ (I-1) *W 
TO CONTINUE 


DO 15 I=1,NUMY 
Y(I)=(I-1)*SPACE+(I-1)*L 
το CONTINUE 


NUMEL=NUMX* NUMY 
WRITE (5,98) NUMEL 
98 FORMAT (I3) 


K=0 
N=0 
ORIENT=0 
me 2ο I=1,9 
DO 22 Jz1,4 

K=K+1 
WRITE(5,99)X(J),Y(I),SMAGD(I,J),SPHASED(I,J),N 


22 CONTINUE 
20 CONTINUE 


99 BSEMATUOES.3,1X,F8.3,1X,F8.3,1X,F8.3,1X,I1) 


STOP 
END 
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APPENDIX G 
I. RADIATION PATTERNS 
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Figure G-1. Beam 1 Radiation pattern (=90°) 
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Figure G-2. Beam 2 Radiation Pattern (ΦΞ905) 
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Figure G-3. Beam 3 Radiation Pattern ($290?) 
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Figure G-4. Beam 4 Radiation Pattern ($290?) 
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Figure G-5. Beam 5 Radiation Pattern (ΦΞ905) 


(dB) 


E 
“4 
© 
O 
oO 
> 
-e 
42 
© 
~ 
ο 
m. 


-75 -60 -45 30 15 109805525011 09259062 08 PES 
Theta (deg) 





Figure G-6. Beam 6 Radiation Pattern (04-90?) 
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II. PATTERN SEPARATION 
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Figure G-7. Simultaneous Radiation 
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Figure G-8. Simultaneous Radiation 
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Figure G-9. Simultaneous Radiation 


III. ANTENNA FOOTPRINT 
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Figure G-10. Beam 1 Relative dB Gain Footprint 
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Relative dB Contours 
(Beam scanned 27 degrees from nadir) 
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Figure G-11. Beam 2 Relative dB Gain Footprint 
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Figure G-12. Beam 3 Relative dB Gain Footprint 
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Relative dB Contours 
(Beam scanned 10 degrees right of nadir) 
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Figure G-13. Beam 4 Relative dB Gain Footprint 
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Figure G-14. Beam 5 Relative dB Gain Footprint 
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Relative dB Contours 
(Beam scanned 45 degrees from nadir) 


Forward 


κα 
O 
© 
H 
[1 
B 
ΓΗ 


ο ος 239 το. ο 15 30 9459.60 75 90 


Left Right 
Cross Track 





Figure G-15. Beam 6 Relative dB Gain Footprint 
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